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We present total cross sections for single and double ionization of helium by antiproton impact 
over a wide range of impact energies from 10 keV/amu to 1 MeV/amu. A non-perturbative time- 
dependent close-coupling method (TDCC) is applied to fully treat the correlated dynamics of the 
ionized electrons. Excellent agreement is obtained between our calculations and experimental mea- 
surements of total single and double ionization cross sections at high impact energies, whereas for 
lower impact energies, some discrepancies with experiment are found. At an impact energy of 1 
MeV we also find that the double-to-single ionization ratio is twice as large for antiproton impact 
as for proton impact, confirming a long-standing unexpected experimental measurement. 
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The double ionization of helium by ion impact has 
been a long and fruitful field of study in atomic col- 
lision physics. As a fundamental four-body system, it 
provides stringent tests of any theoretical description of 
charged particles moving in a Coulomb field. Moreover, 
the sign and magnitude of the ion charge can be exper- 
imentally varied relatively simply to explore interesting 
physical effects. New antimatter collision experiments 
are planned at FAIR, the Facility for Antiproton and Ion 
Research 0. FAIR is an international collaboration on 
atomic and molecular physics that intends to investigate 
antiproton driven ionization processes and even kinemat- 
ically complete antiproton collision experiments. These 
experimental efforts complement the recent intense activ- 
ity in antihydrogen studies of recent years 0] . In the last 
20 years, experimental measurements showed, unexpect- 
edly, that the ratio of double-to-single ionization of he- 
lium from proton impact was around a factor of two lower 
than that from antiproton impact a feature which 
has not yet been observed in a fully converged theoreti- 
cal calculation. Subsequent experimental measurements 
Q of the total double ionization cross section by antipro- 
ton impact not only revealed a larger cross section than 
for proton impact, but also that at low impact energies, 
the double ionization cross section for antiproton impact 
does not monotonically decrease, but shows an increasing 
cross section as the impact energy is lowered. 

In response to this renewed experimental activity, in 
this Letter we present fully converged calculations of sin- 
gle and double ionization of helium from antiproton im- 
pact and find similar unexpected ratios of the double-to- 
single ionization of helium compared with proton impact. 
We also provide converged single and double ionization 
cross sections for antiproton impact over a wide energy 
range from 10 keV to 1 MeV. Theoretical calculations 
for single ionization of helium by antiproton impact cross 
sections have suggested that the experimental measure- 



ments [5] below 30 keV may be inaccurate 0, Q. Ex- 
perimental efforts are underway to remeasure the single 
and double ionization cross sections 0] and preliminary 
results [8| suggest that the single ionization cross section 
does not decrease as rapidly as reported by Hvclplund et 
al These preliminary results may bring closer agree- 
ment between theory and experiment for single ioniza- 
tion. 

For double ionization of helium by antiproton impact, 
no theoretical description has been used to calculate the 
cross section over the full energy range of the measure- 
ments. Barna et al [9i] calculated the double ionization 
cross section at a single energy of 3.6 MeV which agreed 
well with experiment. A close-coupling method was used 
by Di'az et al ^lOj to calculate double ionization cross 
sections for impact energies above 200 keV, with limited 
success. The calculations reported were around a factor 
of two lower than the measurements of [4] . Keim et al cal- 
culated the single and double ionization of helium by an- 
tiproton impact using the framework of time-dependent 
density- functional theory [ll|. These results produced 
reasonable agreement for the single ionization cross sec- 
tions. However, the double ionization cross sections were 
considerably higher than the experimental results for all 
energies. A similar approach was used by Tong et al jT3 |. 
who calculated single ionization cross sections that were 
lower than the experimental measurements, and double 
ionization probabilities for 15 keV and 100 keV that yield 
much larger cross sections compared with experiment. 
The multicut forced impulse method (mFIM) has been 
the only calculation for double ionization of helium by an- 
tiproton impact below 200 keV [l3| which has compared 
well with experiment. However, at these low impact en- 
ergies the mFIM method predicts single ionization cross 
sections that are in less than satisfactory agreement with 
previous calculations [1, 0] . 

In this work, we apply our TDCC method to treat 
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the double ionization of helium by antiproton impact 
[T^ . The six-dimensional time-dependent Schrodinger 
equation for the outgoing electrons is reduced to a cou- 
pled set of two-dimensional partial differential equations 
through a partial-wave expansion of the two-electron 
wavefunction in spherical coordinates centered on the 
target atom. The two-electron wavefunction is subjected 
to a time-dependent projectile interaction which is in- 
cluded through a multipole expansion. We assume a 
straight-line motion for the impacting ion. 

Our previous study using this approach explored 
alpha-particle collisions with helium. For high impact 
energies, excellent agreement between theory and exper- 
iment was found for single and double ionization of he- 
lium. At low energies, charge transfer to the impact- 
ing ion may become significant, so that a description us- 
ing a two-electron wavefunction centered on the target is 
unsuitable. However, in the case of antiproton impact, 
charge transfer cannot occur, so our method is appropri- 
ate for treatment of the antiproton- atom collision at any 
impact energy. 

The fully correlated wavefunction, "J/^^^ , for the single 
and double ionization of a two-electron target atom by 
antiproton collision is obtained by the evolution of the 
time-dependent Schrodinger equation in real time: 



dt 
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where the non-relativistic Hamiltonian is given by 
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and Zp is the projectile atomic number and Zt is the tar- 
get atomic number. For straight-line motion, the mag- 
nitude of the time-dependent projectile position is given 
by 



Rit) = \/&^ + {do + vty 
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where h is an impact parameter, do is a starting distance 
{do < 0), and v is the projectile speed. If we expand 
yj,LA/ jj-^ coupled spherical harmonics and substitute into 
Eq. (1), the resulting close-coupled equations for the ra- 
dial expansion functions, P;^*^ (ri , r2 , t) , are then given 
by 
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where VJ^^^ (ri,r2) is the electron-electron potential 

and W^f^'i^^,^^ {ri,R{t)) represents the interaction be- 
tween the projectile and the electrons. The initial value 
boundary condition for Eq. (4) is given by 

PtJ'{ri,r2,t = 0) = (5L,Lo<5A./,MoP^?f °(^i'^2,T ^ oo) , 

(5) 

where P^^i^^" is the radial portion of the ground state 
wavefunction. The two-electron ground state wavefunc- 
tion of the helium atom is obtained by relaxation of the 
time-dependent Schrodinger equation in imaginary time 
(r —it). The fully correlated ground state is expanded 
in coupled spherical harmonics resulting in a set of close- 
coupled equations for the -P;^°*^° (''i , f 2 , t) radial expan- 
sion functions with the initial boundary conditions given 
by a product of single particle bound radial orbitals for 
the one-electron target ion with L ~ M = . 

We solve the TDCC equations, Eq. (4), using a 



lattice technique to obtain a discrete representation of 
the radial expansion function and all operators on a 
two-dimensional grid. The single and double ioniza- 
tion probabilities, 'P(£', b), are calculated using the time- 
dependent radial wavefunction propagated for a suitable 
time. The radial wavefunction is then projected onto 
both bound and continuum single particle orbitals, as 
discussed in detail in . 

The TDCC calculations employed a 384 x 384 point 
radial lattice with a uniform mesh spacing of Ar = 0.20, 
giving a converged ground state of helium on the 
lattice. In order to achieve accurate double ionization 
cross sections it was found that up to 101 coupled 
channels had to be included in the expansion of the final 
state two-electron wavefunction. The time-dependent 
projectile potential was expanded into monopole, dipole, 
quadrupole, and octopole terms, all of which were 
required to converge the total cross sections. The 
time-dependent wavefunction was propagated from a 
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time when the projectile was a distance do = — 50ao, 
through closest approach, to a time when the projectile 
was at a distance do = +40ao. The total cross sections 
calculated at the lowest impact energy of 10 keV had 
to be propagated for 142a. u. of time to reach this dis- 
tance using a number of time steps greater than 28,000. 
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FIG. 2: Evolution of the ionization probability, 'P(i?, b — 
0.5ao), for a 50 keV antiproton collision with a helium atom 
as a function of the impacting ion distance. Upper panel; 
single ionization probability summed over all partial waves; 
lower panel: double ionization probability. The helium atom 
is located at the origin of the collision system {do — 0). 



FIG. 1: (a) Cross sections for single ionization of helium by 
antiproton impact. Filled squares: TDCC calculations; filled 
diamonds: calculations of Schultz and Krstic Ja; filled and 
open circles: experimental measurements of (b) Cross 

sections for double ionization of helium by antiproton impact. 
Filled squares: TDCC calculations; crosses: calculations of 
Diaz et aZ[l3|; filled circles: experimental measurements of 
i- 



In Fig. 1, we present our TDCC calculations for sin- 
gle and double ionization of helium by antiproton colli- 
sions. The top panel of Fig. 1 shows the single ionization 
cross section measurements [1, HBl, our TDCC results 
(red squares) , and the calculations of Schultz and Krstic 
0. The two sets of calculations are in good agreement 
over a wide energy range down to around 30 keV. Be- 
low 30 keV, the TDCC calculations are higher than the 
experimental measurements, although new preliminary 
measurements of this single ionization cross section [^] 
suggest that the previous measurements 3] may be too 
low. 

The lower panel of Fig. 1 shows the double ionization 
of helium by antiproton impact for ion energies from 10 
keV to 1 MeV. The TDCC (red squares) calculations 
are in excellent agreement with the experimental mea- 
surements above 20 keV. Between 10 and 20 keV, the 
TDCC calculations decrease monotonically. However, 
the experimental measurements show a sharp increase 
as the impact energy decreases. The TDCC results 
are in excellent agreement with the experimental data 
points at 15 keV and 20 keV, suggesting that per- 
haps only the lowest energy measurements are in error. 

A time-dependent propagation of the two-electron 
wavefunction can be a powerful tool in understanding 
the evolution of the collision system. Fig. 2 shows the 
ionization probability, Vi^E, b) for both single and double 



ionization for a 50 keV antiproton collision with a helium 
atom. The initial ion starts at do = — 50ao, then prop- 
agates through distance of closest approach (6 = O.Soo) 
to a final distance of do — -l-40ao. The total propagation 
time for the 50 keV collision system is 63a. u. of time. 
The helium atom is placed at do = Ooo. Fig. 2 shows 
the evolution of the ionization probability for one im- 
pact parameter, {b = 0.5ao) and for an ion distance from 
do = — 5ao — > 15ao. The choice of the impact parame- 
ter equal to 6 = 0.5ao represents the impact distance of 
maximum double ionization probability. The total single 
and double ionization cross sections for a given energy 
are obtained by the relation 

/>oo 

a{E) = 27T P{E, b)bdb . (6) 

The ionization probabilities for both single and double 
ionization of helium by antiproton impact increase 
rapidly just before the antiproton reaches the helium 
atom (positioned at do = Oao), and then tend to a 
constant value beyond do ~ lOao. The results in 
Fig. 2 show that the antiproton collision with the 
helium atom occurs in approximately tcoiusion = 4.24 
a.u. of time or 0.1 fs. Such fast collision times are 
currently just beyond the pulse lengths of experimental 
ultrafast laser techniques which are being considered for 
use in probing electronic behavior of atomic systems. 

The total ionization cross sections are calculated by 
integrating the ionization probabilities as defined in 
Eq. (6). Fig. 3 shows the weighted probabilities (the 
integrand of Eq. (6)) for both proton (black dashed line) 
and antiproton (red solid line) ionization as a function 
of the impact distance from the helium atom at an in- 
cident energy of 1 MeV. For the antiproton impact, the 
peak of the single ionization weighted probability is at 
an impact parameter oi b k, l.Oao, while the double ion- 
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FIG. 3: Weighted probabilities for single (upper) and dou- 
ble (lower) ionization of helium by antiproton (red line) and 
proton (black dashed line) impact at 1 MeV impact energy. 



ization weighted probability peaks at a smaller impact 
parameter (6 « O.Sag). The proton's maximum ioniza- 
tion probability is similar to the antiproton for single 
ionization but shifted to smaller impact parameters for 
double ionization. However, of much greater significance 
is the difference between the proton and antiproton dou- 
ble ionization probabilities. As also observed in the ex- 
perimentally measured double-to-single ratios , our 
calculations find that for antiprotons the double ioniza- 
tion cross section is larger than the proton double ioniza- 
tion by approximately a factor of 2, but that the single 
ionization processes have approximately the same total 
cross section. At this impact energy, the total single 
and double ionization cross sections are in good agree- 
ment with experiment for antiproton impact (as shown in 
Fig. 1) and for proton impact (compared with the mea- 
surements of The reasons for this difference be- 
tween the double ionization cross sections is still not well 
understood. Clearly, the change of sign of the projec- 
tile affects the double ionization mechanism in a much 
greater manner than the single ionization mechanism. 
This could be partially explained by arguments arising 
from a Born-approximation description [l7j] , where inter- 
ference between various mechanisms for double ionization 
depended on the sign of the ion charge, resulting in dif- 
ferent double-to-singic ratios for protons and antiprotons. 
We note that following a rapid ejection of a single elec- 
tron, which based on the upper panel of figure 3 is equally 
probable for the proton impact and antiproton impact, 
the remaining transient collision complexes are quite dif- 



ferent. For proton impact, the remaining electron sees 
a positive charge of -1-3, while for the anti-proton im- 
pact the remaining electron sees a positive charge of -1-1. 
Thus, the second electron may find it easier to escape 
in the antiproton impact case. Similar effects have been 
observed for double photoionization of helium by large 
photon energies, where the shake-off mechanism becomes 
dominant at high photon energy [IS, JLQJ . Although com- 
plicated by the multipole potential is this similar 
mechanism may be responsible for the large double ion- 
ization probability from the antiproton impact. Also, the 
impact distance, 6, can provide useful insight into under- 
standing the dynamics involved in the collision. Clearly, 
double ionization is most probable when the antiproton 
significantly penetrates the electron cloud of the helium 
atom. This suggests the possibility that the antiproton 
undergoes multiple collisions with the electrons to doubly 
ionize. Single ionization is most probable at a somewhat 
larger impact parameter, and still is probable even when 
the ion is three or four atomic units from the atom, at 
which distance there is almost no probability of double 
ionization. Thus the single ionization mechanism is likely 
dominated by a more standard 'binary' collision process. 

In conclusion, the TDCC method has been shown to 
accurately compute single and double ionization total 
cross sections for antiproton impact ionization of helium 
over a wide energy range. The single ionization total 
cross section calculations are in excellent agreement with 
experiment and with previous calculations. The double 
ionization total cross section calculations also agree well 
with experiment above 20 keV. The antiproton double 
ionization cross sections are 2 times larger than the pro- 
ton double ionization cross sections at the same impact 
energy whereas, the single ionization cross sections are 
similar. We aim to further exploit our method to ex- 
tract the differential cross sections for the two outgoing 
electrons after ion impact. Such calculations may allow 
comparison with the large set of experimental measure- 
ments of differential cross sections for ion impact which 
currently exist (e.g. [13]), and should yield more insight 
into the nature of the double ionization process. 

A portion of this work was performed under the 
auspices of the US Department of Energy through Los 
Alamos National Laboratory and through DOE and NSF 
grants to Auburn University. Computational work was 
carried out at the National Center for Computational Sci- 
ences in Oak Ridge, TN, and using the Institutional Com- 
puting Resources at Los Alamos National Laboratory. 
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